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Abstract 

We review a model-based rather than phenomenological approach to low- 
temperature anomalies in glasses. Specifically, we present a solvable model inspired by 
spin-glass theory that exhibits both, a glassy low-temperature phase, and a collection 
of double- and single-well configurations in its potential energy landscape. The distri- 
bution of parameters characterizing the local potential energy configurations can be 
computed, and is found to differ from those assumed in the standard tunneling model 
and its variants. Still, low temperature anomalies characteristic of amorphous mate- 
rials are reproduced. More importantly perhaps, we obtain a clue to the universality 
issue. That is, we are able to distinguish between properties which can be expected 
to be universal and those which cannot. Our theory also predicts the existence, under 
suitable circumstances of amorphous phases without low-energy tunneling excitations. 



1 Introduction 

Ever since the first measurements of Zeller and Pohl ffl revealed that the specific heat 
and and the thermal conductivity of glassy systems at low temperatures are strikingly 
different from those of crystalline substances, the low temperature physics of glassy and 
amorphous materials has been a subject of intense research efforts, both experimentally 
and theoretically. 

Specifically, it was found U that below approximately 1 K the specific heat of glassy 
materials scales approximately linearly with temperature, C ~ T while the corresponding 
scaling for the thermal conductivity k is approximatly quadratic k ~ T 2 . Both findings 
contrast the T 3 behaviour of these quantities in crystals. Between 1 K and approximately 
20 K the thermal conductivity exhibits a plateau and, then continues to rise as the te- 
perature is further increased. The specific heat, too, changes its behaviour in the 1-20 K 
regime. It exhibits a bump if displayed in C/T 3 plots. 

*To appear in: Festkorperprobleme/ 'Advances in Solid State Physics 38 (1998) 



1 



Particularly intriguing is the fact that the anomalous behaviour below 1 K appears to 
be universal in the sense that it is shared by virtually all glassy and amorphous materials, 
whereas the behaviour between 1 and 20 K seems to display a greater material dependence. 
The universality in the thermal conductivity data actually appears to go beyond the level 
of exponents, in that by a suitable rescaling with a constant depending on the Debeye 
temperature and the the sound velocity, thermal conductivity data for various substances 
can be scaled onto one master-curve both below and above, though not on the plateau 
I- 

It should be added that there is also an extensive body of experimental data char- 
acterizing the anomalous response of glassy and amorphous materials to external probes 
such as ultrasound or electromagnetic fields. The reader is invited to consult the reviews 
of Hunklinger et al. || or Phillips [|| on these matters. 

The unusual material properties at low temperatures are usually attributed to the 
existence of a broad range of localized low-energy excitations in amorphous systems - 
excitations not available in crystalline materials. At energies below 1 K, these are believed 
to be tunneling excitations of single particles or (small) groups of particles in double- 
well configurations of the potential energy (DWPs). This is the main ingredient of the 
phenomenological so-called standard tunneling model (STM), independently proposed by 
Phillips || and by Anderson et al. As a second ingredient of the STM, it is supposed 
that the local DWPs in amorphous systems are random, and specific assumptions concern- 
ing the distribution of the parameters characterising them must be advanced (see below) 
to describe the experimental data below IK || ||. 

Neither Phillips || nor Anderson et al. Q consider excitations other than the tun- 
neling excitation in the DWP structures postulated by them — their local degrees of 
freedom are thus true two-level systems (TLSs). As a consequence their model is unable 
to account for the different physics that appears in the temperature range between 1 and 
20 K. However, different sets of low energy excitations in amorphous systems can exist: (i) 
Higher excitations in the abovementioned DWPs — a possibility that is surely contained 
in the physical picture advanced by Phillips or Anderson et al., but was justly consid- 
ered irrelevant for the very low temperature regime and then somehow not reconsidered 
when problems with the model arose at temperatures above 1 K. (ii) Localized vibrations 
(localized phonons) in soft anharmonic single well configurations of the potential energy. 
The latter have been postulated (in addition to tunneling excitations in DWPs) within 
the, likewise phenomenological soft-potential model (SPM) 0, H to describe the physics 
above 1 K. Again, local potential energy configurations are supposed to be random, and 
specific assumptions concerning the distributions of the parameters characterizing them 
are required to account for the experimental data also above 1 K, such as the crossover to 
T 3 -behaviour of the specific heat and the plateau in the thermal conductivity Jj], B|. 



Both, localized soft vibrations || and DWPs [10| responsible for two-level tunneling 
systems (TLS) have been seen in molecular dynamics studies of Lennard Jones glasses. In 
the case of soft vibrations || no attempt, however, was made to determine the shape of the 
single-well potentials (SWPs) supporting these vibrations, so as to check the hypotheses 
of the SPM. On the other hand, local potential energy configurations giving rise to TLSs 
were analysed within the confines of a generalized SPM, which assumes that locally 
the potential energy surface (along some reaction coordinate) can be described by certain 



2 



fourth order polynomials, with coefficients distributed in a specific way. No single-well 
configurations were taken into account, though, to determine the statistics of the coeffi- 
cients |10|| , as in principle they should to make full contact with the assumptions of the 
SPM. 

Some time ago, in a lucid critical discussion of the standard theory of TLSs in amor- 



phous systems, Yu and Leggett [11] pointed out that there was nothing in the STM that 
could reasonably account for the considerable degree of universality observed in amor- 
phous systems. Analogous remarks would mutatis mutandis apply to the SPM, and it is, 
of course, mainly related to the phenomenological nature of these models. They have to 
rely on assumptions concerning in particular the distributions of parameters character- 
ising local potential energy configurations, which — while plausible in some respects - 
are certainly much less so in others, and are lacking support based on more microscopic 



approaches, such as that of [12] for KBr:KCN mixed crystals. Neither model, to be sure, 
accounts for a mechanism that would explain how the required local potential energy 
configurations would arise, and how they would do so with the required statistics. 

It is here that our model-based approach to glassy low-temperature physics [13, 14], 
which was started about two years ago, attempts to fill a gap. 

Our approach takes as its starting point the very observation of universality of glassy 
low-temperature anomalies, which we translate into a strategy as follows. Since — in the 
light of universality — the detailed form of particle-interactions is apparently to a large 
extent immaterial to the phenomena that are to be modeled, we should be justified in 
proposing something like a caricature-glass. We would like this to be understood here in the 
sense of something as simple as possible as long as it retains essential properties of glass- 
forming systems. Our success will of course depend on the quality of our understanding 
as to what these essential ingredients might possibly be. We take them to be (i) particles 
moving in continuous space, subject to (ii) an interaction (any interaction!) that gives 
rise to an amorphous low temperature phase. Within these confines, our choice of the 
interaction has been mainly guided by the demand for simplicity, and analytic tractability. 
The outcome of these deliberations has been a proposal inspired by spin-glass theory 



1 13, 14] that describes a glassy material as an anharmonic interacting particle system with 
random interactions at the harmonic level, the details of which will be specified in Sect. ^ 
below. The interactions are chosen in such a way that the model is solvable via mean-field 
and replica techniques well konwn in the theory of spin-glasses [ O ] . 

We have organized the remainder of our material as follows. In Sect. ||] we present 
details of our model, and review main features of its solution within mean-field theory. 
We compute its glass transition temperature and exhibit its phase diagram, featuring 
ergodic, polarised and glassy phases. Sect. || is devoted to mapping out the potential energy 
surface of the system in (one of) its classical glassy ground-states, and to determining the 
statistics of its local potential energy configurations. Thermodynamic consequences at 
low temperatures are explored in Sect. El where we look at excitation spectra in local 
potential energy configurations, at the density of states, and at the specific heat. Dynamic 
consequences are briefly considered in Sect. ||. We close in Sect. |6| with a summary and 
with an outlook on open issues. 
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2 Spin— Glass Approach and a Solvable Model 



We suggest to consider the following Hamiltonian as a candidate for the description of 
glassy low-tmeperature properties 

N 2 

W = E|^ + C/ int(M) , (1) 

r— f 2m 

i=i 

with an interaction energy given by 



hit 



({«<}) = - J E + - E G 0>o . ( 2 ) 



in which we include an on-site potential of the form 

ml \ a 2 2 , a 4 4 fo\ 

G(v) = —v +—v . (3) 

For the time being we specialize to 02 = 04 = 1. 

The description is in terms of localised degrees of freedom, i.e., the V{ may be interpreted 
as deviations of particle positions from a given set of reference positions — in a spirit akin 
to that used to describe dynamical properties of crystalline solids. Thus, we assume the 
system to be already in a solid state, and do not attempt to provide a faithful description 
of the liquid phase.Q 

We propose to model the glassy aspects by taking the interaction at the harmonic level, 
that is, the JJy, to be random, so that the reference positions will generally be unstable 
at the harmonic level. This is why we have added the on-site potentials G(v), namely 
to stabilize the system as a whole. By including in G a harmonic term, we can use the 
parameter 7 to tune the number of modes that actually are unstable at the harmonic level 
of description. Apart from its role in this 'tuning' aspect, we have opted for the simplest 
choice that achieves the required stabilization. 

We believe our opting for random interactions to be justified by appeal to universality 
(the detailed form of the interactions cannot be crucial for the emergence of glassy low- 
temperature anomalies), and by the observation made in recent years of the fundmental 
similarity between truly quenched disorder and so-called self-induced quenched disorder 



We choose the random interactions in such a way that the system is amenable to 
analysis within replica mean-field theory, well known from the theory of spin— glasses fl5|| . 
Specifically, we choose them to be Gaussian random variables with Jy = Jq/N and Jf — 
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Jij = 1/N. The harmonic part of J7int({^i}) is thus reminiscent of the SK spin-glass [17], 
except for the fact that we are dealing here with continouus degrees of freedoms rather 
than with classical Ising spins. Once more, by appeal to universality, and by noting that 



1 Note that we are writing down what appears to be a scalar model here. The 3-D nature of particle 
coordinates is, however, included in our description, if we choose to consecutively enumerate the cartesian 
components of the particles, in which case N has to be read as three times the nuber of particles of the 
system. 
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we do not intend to describe correlations at a critical point, it is not unreasonable to argue 
that a mean-field description should be sufficient to reveal the essential physics we are 
after. 

Let us mention for completeness that models with an interaction energy of the form (|2|) 
were considered before in an entirely different field, namely in the context of analog-neuron 
systems and networks of operational amplifiers [jl8| , |l9| . Incidentally, this connection opens 
up an unexpected but rather fascinating perspective of studying glassy dynamics by em- 
ulating it directly in highly integrated circuitry, albeit perhaps not in the regime of low 
temperatures, where quantum effects are important. 

At the classical level of description, all interesting information about the thermody- 
namics of the system is obtained from the configurational free energy 

f N (P) = -(PN^ln [ Y[dvieM-mnt({vi})} . (4) 
J i 

The free energy has to be averaged over the ensemble of possible relisations of the Jij 
configurations so as to get typical results. This is accomplished by means of the so-called 
replica method, well known in the theory of disorderd systems (see, e.g. |]15|). We will not 
here reproduce the calculations, but merely quote the result. The quenched free energy 
- that is, the disorder-average of (||) — is represented as the limit f((3) = lim n ^o fn(P), 
with 

nfn(P) = i/o£^ + |/?£<4 

a. a,b 

- /T'ln / n^ a exp[-/3C/ cfr (K})] • (5) 

Here 

U eS = - Jo E Pa v a - ^ £ q ab v a v b + G(v a ) (6) 

a a,b ^ a 

is a replicated single-site potential and the order parameters p a = N~ x J2i i v t) an d Qab 
X 1 J^i ( v f v i) are determined from the fixed point equations 

p a = (v a ) , a = l,...,n (7) 
q a b = (v a v b ) , a,b = l,...,n , (8) 

where (...) denotes a Gibbs average corresponding to the replica potential (g). The limit 
n — > is eventually to be taken in these equations. 

Eqs. (P)-(0) were analysed in the replica symmetric (RS) and the 1st step replica- 
symmetry breaking (1RSB) approximations fl4|| . In RS one assumes p a = p for the 'po- 
larization', and q aa = q and q a b = q for a ^ b for the entries of the Edwards- Anderson 
matrix. These are given as solutions of 

P =((v)h, C ={A=^-) Z > v = ((v) 2 h- (9) 
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Figure 1: (a): Glass temperature T g as a function of 7 for Jo = 0. (b): T = phase 
diagram in RS. E denotes the T = 0-limit of the ergodic phase, G is the glassy phase, 
and P is a phase with macroscopic polarization. The full line is the AT line (see text). 
Below the dot-dashed line is the region with DWPs, as described in detail in Sect. [|. 
The short-dashed line separates the glassy phase G from the phase P with macroscopic 
polarization. From fl4| 



Here C = j3{q — q), and (. . .) z denotes an average over a standard Gaussian z while (. . .) 
is now a thermal average corresponding to the RS single-site potential 

U eS (v) = -h eS v - X -Cv 2 + -G(v) , (10) 
2 7 

with 

Kb = hns = JoP + y/qz (11) 

It turns out that our model exhibits a transition from an ergodic phase with p = q = 
to a glassy phase with q / at some temperature T g which depends on Jo and 7. For 
sufficiently large Jo, the transition is into a polarized phase withp 7^ 0. The glass-transition 
temperature as a function of 7 is shown in Fig. [l|a. In Fig. [l|b, we exhibit also the T = 
limit of the phase diagram, which will become of importance for the discussion of of the 
potential energy surface of the model in Sect. || below. 

The assumption of RS is not justified at low temperatures and large 7, where replica 
symmetry breaking (RSB) — an instability that breaks the permutation invariance among 
the n replica — is expected to occur. The location of the instability against RSB is given 
by the de Almeida-Thouless (AT) criterion @ 1 = (3 2 (((v 2 ) -(v) 2 ) 2 ) z . In a one- 
step replica-symmetry breaking (1RSB) approximation, which is believed to constitute 
a major step towards the full solution, U e s(v) is of the same form (|To|), however with 
h RS replaced by h 1RSB = J p + z + v% - 9o^l, and C by C = (3{q - q x ), where 
we use standard notation [15, 17] for the entries of the (/-matrix. Along with a so-called 



partitioning parameter m, they are determined from a more complicated set of fixed point 



equations ]15| , 21| . In 1RSB, zq is a Gaussian, whereas the distribution of Z\ is more 
complicated. 
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In the following section we will set out to demonstrate that the mean-field solution of 
the model we have obtained here does in fact also contain information about the potential 
energy surface of the model, the statistics of which is generally believed to be crucial for 
the emergence of glassy low-temperature anomalies. 

3 Mapping Out the Potential Energy Surface 

To map out the potential energy surface of the system, we observe that all that mean field 
theory is about, is to represent the interaction energy of the system as a sum of effictive 
independent single site potential energies 

U int ({vi}) — ^C/efffa) , (12) 

i 

self-consistently to be determined in such a way as to get thermal averages correct. In 
the context of non-random systems, the most famous example is the Curie- Weiss model of 
ferromagnetism, where the effective single site potential is the Zeemann energy of a spin 
in an effective (non- fluctuating) local- or mean-field determined by its neighbours. 

In our case, we obtain the sum ( |i~2| ) of independent single site potentials U e g (vi), which 
contain random parameters. Replica theory achieves nothing but a selfconsistent deter- 
mination of the distribution of these random parameters. Within the RS approximation 
these effective local single site potentials are given by (|IT}),(|rI]), containing a single random 
parameter, viz the Gaussian distributed effective fields /jrs, having mean Jqp and variance 
q. The parameters of p, q and C characterizing the Gaussian ensemble of single-site po- 
tentials U e ff(v) are determined self-consistently through (P)-([TT|). In 1RSB, the ensemble 
of single site potentials is again of the form fllPP, but now with /irs replaced by the more 
complicated /iirsb and C by C = f3(q — q%) as noted above. Again there is only a single 
random parameter, a locally varying effective random field. This feature is easily seen to 
persist at all levels of replica symmetry breaking. 

There is one additional ingredient, namely we take the the T = limit of the theory 
to select one of the (possibly many) collective classical glassy groundstate configurations. 

In terms of the mean field solution, we thus have a representation of the glassy potential 
energy landscape as an ensemble of randomly varying effective single-site potentials. Their 
statistics is determined by the statistics of the effective fields /irs and ^-irsb in the RS 
and the 1RSB approximation, respectively. 

Now note in particular, that the initial stabilising concave upward single site potential 
G(v) also appears in U e g(v), but there is now an additional harmonic term, viz. — \Cv 2 
- entirely of collective origin — that renormalises the effective total harmonic restoring 
force. For 

C-->0 (13) 

7 

the total harmonic contribution to U e g (v) becomes convex downward near the origin v = 0, 
so that for sufficiently small /irs/irsb the effective single site potential U e g(v) attains a 
DWP-form, which is of collective origin! The region in parameter space where DWPs exist 
and its boundary are indicated in the T = phase diagram in Fig. |l]b. 
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We are now able to make contact with the phenomenology of glassy low-temperature 
anomalies, and to compare with the ideas underlying their explanation within the STM 

|3| or the SPM [jTj. The main difference and, we believe, important new perspective 
added by our approach is that the appearance of DWPs and the statistics of local potential 
energy configuraqtions need not be hypothesised from the outset, but rather that it emerges 
as a collective effect, originating in the frustrated interactions Jj,-, if the temperature is 
sufficiently low and the parameter 7 sufficiently large. Moreover, the statistics of the local 
potential energy configurations has been rendered part of the world of the computable. It 
follows directly from the statistics of the /irs/irsb- 

To explore the consequences for the low temperature behaviour of glasses quantita- 
tively, one solves the quantum mechanical problem for U e fi(vi) described be the Hamilto- 
nian 

#eff = |~ + U eS (v) (14) 

with U e g(v) given by (|io|) , (|ll|) in the RS approximation, or by the corresponding expression 
valid for the 1RSB (or higher order RSB) ansatze. 

In the STM one approximates the tunnel-splitting e of the ground-state in the DWP 
case to be given in terms of the asymmetry A of the DWP and the tunneling-amplitude 
A by the TLS-result 



A 2 + Ag , (15) 

the value of Ao being related with the barrier-height V between adjacent minima and 
the distance d between them. In WKB-approximation one has Ao — TujJq exp(— A), with 
A = f {2mV/h 2 ) 1 / 2 , where uq is a characteristic frequency (of the order of the frequency 
of harmonic oscillations in the two wells forming the double well structure), and m the 
effective mass of the tunneling particle. Moreover, a specific assumption is advanced con- 
cerning the distribution of asymmetries A, and tunneling-matrix elements Ao, namely, 
in terms of A and A one assumes -P(A, A) ~ const. Note that this assumption implies in 
particular also that these quantities are uncorrelated. 



We can check these assumtions against the results from our model [13]. Obviously, since 
the values of both A and A are determined by a single random variable, viz. /irs/irsbj these 
quantities must be perfectly correlated in our model, in contrast to the assumptions of the 
STM, although their distributions separately are indeed rather flat, except for anomalies 
at their upper end (see Fig. |2|) where due to large ^Irs/irsb DWPs cease to exist. 

The perfect correlation can be weakended, but not eliminated, by allowing local ran- 
domness, that is, by making the coefficients 02 and 04 in (||) vary randomly from site to 
site. A randomly varying a% is, in fact, a natural choice in our approach, since the a% con- 
tribution to ZJjtnt should be considered part of the harmonic level, which we had assumed 
to be random to begin with. Indeed, by allowing this kind of randomness, our results can 
be brought into better agreement with the experimental situation (see below). As to 04, 
we have observed that taking it to be random (within limits of course) has no dramatic 
effect, so we adhere to it remaining non-random, serving solely, as it should, its stabilising 
purpose. We should add that our model remains solvable with these modifications. Outer 
averages in the fixed point equations simply have to be read as implying an additional 
average over the local randomness. 
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Figure 2: Distributions 'P(A) (for positive A), and 'P(A) in RS. Here 7 = 4, while Jo = 0. 
After lO| 



As a final unexpected result, note that our theory predicts the existence under suit- 
able circumstances of an amorphous phase without DWPs, and hence without low-energy 
tunneling excitations (see Fig. |l]b). Until very recently, this would perhaps have univer- 
sally been considered rather a surprise. However, recent internal friction measurements in 



specially prepared amorphous Si:H films by Liu et al. |22|] indicate that such a possibility 
has to be seriously taken into account. 



4 Thermodynamics 



Information about the thermodynamics of the system at low temperatures, where it is 
dominated by quantum effects, is obtained by studying the excitation spectra of the local 
Hamiltonians (14), and by using this information — avaraging it over the randomness char- 
acterizing the ensemble of these local Hamiltonians, i.e. over the /iRg/iRSB _ distributions 
(and possibly the (^-distribution) — to compute densities of state and thermodynamic 
functions the usual way. 

In Fig. ||a we exhibit a typical excitation spectrum, as it varies with the value of h e g. 
The salient features are two kinds of branches, steeper and less steep ones. The former are 
related to tunneling (inter-well) excitations whereas the the latter are of intra-well type. 
That is, the latter would also exist if the wells were not communicating via a barrier of 
finite height and width. Note that these families of levels do not cross as might appear 
at first sight. Rather, there are avoided 'would-be' level-crossings precisely due to the 
tunneling mechanism. At at energies higher than the barrier or at large /irs the avoided 
crossing-pattern disappears, because there is no longer a barrier to be tunneled through. 

The corresponding density of states (DOS) (averaged over the ensemble) is displayed 
in Fig. [3|b. Note its constant value at energies lower than that of the first intra-well 
excitation. It arises because, (i) apart from the immediate vicinity of its lower cutoff Ao, 
the level-splitting of the ground-state scales linearly with the asymmetry A, which in 
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Figure 3: a): Excitation spectrum as a function of /irs- 7 = 4 and Jo = 0. are as in Fig. 
||. b): Density of states in RS, for parameters as in (a). 

turn is proportional to h e $, and (ii) because the distribution of h e g is virtually flat in this 
regime. 

This energy range is to be associated with the universal glassy low-temperature anoma- 
lies below IK! Its properties are mainly driven by the distribution of effective fields h e g 
and are thus entirely of collective origin. In particular it is related to the temperature 
range in which the specific heat scales linearly with temperature (see Fig. 

By contrast, the energies of the intra-well excitations are much less sensitive to the 
value of the effective field, hence less steep in the excitation spectrum if plotted vs. h e ff. 
These give rise to the strong increase of the DOS at the value of the first intra-well 
excitation (the contribution to the DOS being proportional to the inverse slope of the 
level as a function of fo e ff)- The intra-well excitations depend more strongly on the shape 
of the local stabilizing potentials G(v), hence they cannot be of collective origin — basically 
they turn out where we have chosen them to be, with or without local randomness in the 
G(v). So the physics that depends on these features, namely the bump in C(T)/T 3 and 
the plateau in the thermal conductivity, cannot be expected to be universal. 

This separation within our approach of aspects that can be expected to be univer- 
sal from those which cannot is in perfect qualitative agreement with the experimental 
situation. It is not unreasonable to suppose that such a distinction between local and 
global contributions to C/j n t can be made for real materials as well, the contribution of the 
long-range global one being for instance of elastic nature (and well known to give rise to 
frustration). 

Note that our representation of local potential energy configurations has finally turned 
out to be similar to that proposed within the SPM J?], H with some marked differences 
in the details: (i) Within our approach the very occurrence of DWPs is of identifiable 
collective origin, and so is the distribution of asymmetries, (ii) No analogous mechanism 
is available for determining (distribution of) the parameters entering G(v), which are 
mainly responsible for the distribution of barrier heights and the nature of (single- and) 
intra-well excitations which determine the material-specific properties. From our vantage 
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Figure 4: Specific heat at low temperatures a): without local disorder (full line RS, dashed 
line 1RSB) , b): with 7 = 1, and a,2 a Gaussian of mean 1 and variance a = 5 in RS. 
Other parameters are as before. The thin line corresponds to C ~ T . 

point, therefore, one should not even attempt at proposing for these something universally 
valid for all glassy materials. 

There is some influence of local disorder on the behaviour of the low T specific heat also 
in the 'universal' regime changing it to slightly superlinear, with exponents depending on 
the nature of the distribution. This feature, too, agrees nicely with experimental findings, 
and it might be used to give some handle at fixing the distribution of local parameters for 
specific materials. 

Concerning universality, Yu and Leggett supposed that it could result only due to a 
sufficiently long range interaction between TLSs, an idea which has since the been pursued 
in a series of papers by Burin and Kagan ]2^| . Here we have seen that that it can arise 
without assuming interactions at the level of quantized excitations, because the statistics 
of local potential energies is by itself already a largely collective affair. 

5 Dynamics 

Let us add a few brief remarks concerning dynamics. We have in mind here the computation 
of phonon mean free paths, of internal friction and the thermal conductivity, the latter two 
computable once the phonon mean free paths are known, as well as the computation of 
dielectric susceptibilities. All interesting properties can be computed within linear response 
theory. We assume a bilinear interaction of the local coordinates Vi, whose dynamics are 
given by H$ = H e R in (|l4|), with the extended phonon degrees of freedom 




(16) 



ks 



that is mediated by the strain-field, and of the form 




(17) 



s 
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with s labelling the acoustic branches of the phonon spectrum, e s denoting the contribution 
of branch s to the strain field at site i and appropriate coupling constants 7 S . A Debeye 
model is assumed for the phonon bath. 

The dynamic properties of interest can be computed from the imaginary part of the dy- 
namic susceptibility x'vvi u )> which in turn, via fluctuation dissipation theorems, is related 
to the (symmetrized) centered correlation function 

C m (t) = (Av(t)\Av) = -(Av(t)Av + AvAv(t)) , (18) 

specifically to the spectral function associated with it. Here Av = v — (v). The Mori- 
Zwanzig Projection technique and (at this time) a weak coupling assumption are used to 
carry through in practice. 

As of now, we have not yet produced quantitative results. So much can, however, be 
said at this point. In the limit of very low temperatures, our model clearly approaches the 
TLS-physics, and we expect a roughly quadratic temperature dependence of the thermal 
conductivity. The increased density of states in the 1 - 10 K range goes along with a new, 
and faster scattering mechanism: both are due to the existence of intra-well (as opposed to 



the slower inter-well) transitions in this energy range. As discussed by Yu and Leggett [11] 



the combination of these features can be expected to produce the plateu in the thermal 



conductivity. We shall report on these matters in greater detail in the near future [21] 



6 Conclusions and Outlook 

In summary, we have proposed to take a new look at glassy low-temperature anomalies 
from a model-based rather than phenomenological approach. The models so far considered 
do certainly not (yet) attempt to describe the details of any specific substance. Our aim 
has been to formulate something very schematic that should capture the essentials of 
amorphous low temperature phases — a caricature in the same sense as the Ising model is 
a caricature description of uniaxial ferromagnets or, perhaps less boldly, the SK model a 
prototypic description of a spin-glass. Chances are of course, that we still haven't got the 
bare essentials right. But we feel that our first attempts do point in the right direction. In 
particular the clue we have obtained concerning the understanding as to which phenomena 
might be expected to be univeral and which not, does look like an encouraging general 
qualitative result. 

Our approach offers the unique possibility to study the whole range of glassy physics, 
from the regime of the glass transition temperature down to the low temperatures where 
quantum effects play a dominant role, all within a single set of model assumptions. In 
particular, it should be interesting to study the glassy dynamics within the type of models 
we have proposed at or near their respective glass transition temperature. 

Another issue concerns the classification of spin-glasses into roughly two families, 
one exhibiting a continuous transition and infinitely many levels of RSB, the other a 
discontinuous transition, and only one level of RSB, the model class we have presented 
here belonging to the first. There are indications that at least from the point of view of 
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dynamics near the glass transition, the second class should be more appropriate, since 
it exhibits dynamic freezing transitions at temperatures higher than the thermodynamic 
transition temperature seen in equilibrium treatments [24|. We are currently studying a 
modified version of the Bernasconi model |25| in this context [B(| . Finally, we have started 
to apply our approach to the case of defect crystals like KCl:Li to study the effects of the 
dipolar interaction among the Li defects at higher concentration [27], a system that has 
the distinct advantage that the on-site potential G(v) is well known. 

In the present paper we have quantised our system only after a mean-field decoupling. 
In principle, the analysis should start out from a full-fledged quantum-statistical formu- 
lation, using imaginary time path integrals in conjunction with the replica method. Such 
a formulation has been worked out [28| and will be analysed in detail in the near future. 
So much can be said already at this point. Our present approch amounts to investigating 
the quantum-statistical physics in the so-called static approximation |^]. In the present 
context, this approximation amounts to ignoring the feedback of quantum fluctuations on 
the effective single site potentials. 
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